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1.  INTRODUCTION 

This  is  the  final  report  under  Grant  N00014-89-J-1109,  October  1988-Deceniber 
1995.  However,  the  report  also  covers  all  work  done  under  the  predecessor  instrument, 
Contract  N00014-84-K-0574,  August  1984-December  1988.  The  total  period  of  the 
report  is  thus  1  August  1984  to  31  December  1995.  Some  of  the  work  under  Contract 
N00014-84-K-0574  grew  out  of,  or  was  an  extension  of,  projects  from  two  previous 
ONR  contracts,  N00014-75-C-0867  (April  1975-August  1984)  and  N00014-82-K-0805 
(September  1982-November  1984).  Ties  to  the  latter  two  contracts  are  noted. 

The  research  reported  here  was  primarily  in  nonlinear  acoustics.  The  broad  pur¬ 
pose  was  to  investigate  the  behavior  of  high-intensity  sound  waves,  especially  to  study 
phenomena  that  have  no  counterpart  in  small-signal  acoustics.  In  a  number  of  cases, 
however,  research  on  a  topic  in  high-intensity  sound  began  with,  or  was  even  domi¬ 
nated  by,  small-signal  aspects  of  the  problem. 

The  report  is  a  summary.  For  more  details,  see  publications  or  dissertations/theses 
cited  herein  or  listed  in  the  Chronological  Bibliography.  The  following  annual  reports 
give  a  chronological  description  of  the  research: 


•  Contract  N00014-84-K-0574 

1.  1st  Annual  Report,  Period  1,  August  1984-31  October  1985  (85-8)* 

2.  2nd  Annual  Report,  1  November  1985-31  October  1986  (86-6) 

3.  3rd  Annual  Report,  1  November  1986-30  September  1987  (87-4) 

4.  4th  Annual  Report,  1  October  1987-30  September  1988  (88-5) 

•  Grant  N00014-89-J-1109 

1.  1st  Annual  Report,  1  October  1988-30  September  1989  (89-6) 

2.  2nd  Annual  Report,  1  October  1989-30  September  1990  (90-4) 

3.  3rd  Annual  Report,  1  October  1990-30  September  1991  (91-5) 

4.  4th  Annual  Report,  1  October  1991-30  September  1992  (93-2) 

5.  5th  Annual  Report,  1  October  1992-30  September  1993  (93-5) 

6.  6th  Annual  Report,  1  October  1993-30  September  1994  (95-2) 

7.  7th  Annual  Report,  1  October  1994-30  September  1995  (95-3) 

’Numbers  given  in  this  style  refer  to  items  in  the  Chronological  Bibliography  given  at  the  end  of 
this  report,  e.g.,  85-8  means  the  fifth  entry  in  the  list  for  1985. 


1 


This  page  intentionally  left  blank. 


2 


2.  STUDENTS  AND  OTHER  PERSONNEL 


A.  Graduate  Students 

1.  Michael  R.  Bailey 

(a)  Academic  department:  Mechanical  Engineering 

(b)  Student  period:  September  1991-present 

(c)  ONR  support:  September  1991-April  1993,  September  1994-December 
1995 

(d)  Other  support:  Applied  Research  Laboratories  IR&D  Program,  May  1994— 
August  1994 

(e)  Degrees:  M.S.,  May  1994;  Ph.D.,  expected  May  1997 

2.  Charles  E.  Bradley 

(a)  Academic  department:  Mechanical  Engineering 

(b)  Student  period:  September  1987-May  1993 

(c)  ONR  support:  July  1987-June  1993 

(d)  Other  support:  UT  Scholarship,  September  1991-May  1992 

(e)  Degrees:  M.S.,  December  1990;  Ph.D.,  May  1993 

3.  Samuel  C.  Clark 

(a)  Academic  department:  Mechanical  Engineering 

(b)  Student  period:  September  1993-May  1994 

(c)  ONFo  support:  September  1993-January  1994 

(d)  Other  support:  College  of  Engineering  Fellowship 

(e)  Degree:  none 

4.  Frederick  D.  Cotaras 

(a)  Academic  department:  Electrical  and  Computer  Engineering 

(b)  Student  period:  September  1983-December  1988 

(c)  ONR  support:  September  1983-December  1988  (no  salary  support) 

(d)  Other  support:  Canadian  government  (salary  support) 

(e)  Degrees:  M.S.,  August  1985;  Ph.D.,  May  1989 
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5.  Lawrence  J.  Gelin 

(a)  Academic  department:  Mechanical  Engineering 

(b)  Student  period:  September  1993-May  1995 

(c)  ONR  support:  September  1993-May  1995 

(d)  Other  support:  none 

(e)  Degree:  M.S.,  May  1995 

6.  Andrew  J.  Kimbrough 

(a)  Academic  department:  Mechanical  Engineering 

(b)  Student  period:  September  1986-May  1987 

(c)  ONR  support:  June  1986-May  1987 

(d)  Other  support:  NASA 

(e)  Degree:  none 

7.  Ping-Wah  Li 

(a)  Academic  department:  Physics 

(b)  Student  period:  September  1988-August  1993 

(c)  ONR  support:  September  1991-May  1993 

(d)  Other  support:  NIH,  September  1989-August  1991;  Applied  Research  Lab¬ 
oratories  IR&D  Program,  June-September  1993 

(e)  Degree:  Ph.D.,  August  1993 

8.  David  A.  Nelson 

(a)  Academic  department:  Mechanical  Engineering 

(b)  Student  period:  September  1981-December  1984 

(c)  ONR  support:  September  1981-June  1985 

(d)  Other  support:  NASA 

(e)  Degree:  M.S.,  December  1984 

9.  James  A.  Ten  Cate 

(a)  Academic  department:  Mechanical  Engineering 

(b)  Student  period:  September  1985-May  1992 

(c)  ONR  support:  June  1985-May  1992 

(d)  Other  support:  ONR  Grant  N00014-89-J-1103,  Applied  Research  Labo¬ 
ratories  IR&D  Program,  Packard  Foundation,  Texas  Advanced  Research 
Program,  NSF,  Newport  Research  Corp.,  UT  Austin  Univ.  Res.  Inst. 

(e)  Degree:  Ph.D.,  May  1992 
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10.  Youseph  Yazdi 

(a)  Academic  department:  Electrical  and  Computer  Engineering 

(b)  Student  period:  September  1990-present 

(c)  ONR  support:  September  1990-August  1992 

(d)  Other  support:  David  and  Lucille  Packard  Foundation 

(e)  Degree:  none  attributable  to  ONR  support 

B.  Senior  Personnel 

1.  David  T.  Blackstock  —  principal  investigator 

2.  Frederick  D.  Cotaxas  —  After  completing  all  requirements  for  his  doctorate  in 
December  1988,  Cotaras  returned  to  his  permanent  post  at  Defence  Research 
Establishment  Atlantic,  Dartmouth,  NS,  Canada.  He  returned  to  ARL:UT  for 
a  brief  period  in  1991  as  a  consultant. 

3.  Mark  F.  Hamilton  —  After  a  year  at  the  University  of  Bergen  as  the  Acoustical 
Society’s  F.  V.  Hunt  Postdoctoral  Fellow,  Hamilton  rejoined  Nonlinear  Acous¬ 
tics  Division  (NAD)  ARL:UT  in  1984  as  a  postdoctoral  fellow.  In  September 
1985  he  was  appointed  Assistant  Professor  in  the  Mechanical  Engineering  De¬ 
partment,  The  University  of  Texas  at  Austin,  where  he  developed  his  own  ONR 
supported  program  of  research  in  nonlinear  acoustics.  However,  he  continued  to 
collaborate  with  personnel  in  NAD  on  several  of  its  ONR  supported  projects. 

4.  Christopher  L.  Morfey  —  Morfey  (Institute  of  Sound  and  Vibration  Research, 
University  of  Southampton,  UK)  was  an  annual  consultant  during  the  period 
1984-89  and  a  Visiting  Research  Fellow  for  six  months  during  1991. 

5.  Jacqueline  Naze  Tj0tta  —  On  leave  from  the  Mathematics  Institute,  University 
of  Bergen,  Bergen,  Norway,  Naze  Tjptta  was  a  Visiting  Research  Fellow  at 
ARL:UT  during  the  period  1987-90.  Although  only  partly  supported  by  ONR 
funds,  she  contributed  significantly  to  the  projects  of  Cotaras  and  Ten  Cate. 

6.  Sigve  Tj0tta  —  On  leave  from  the  Mathematics  Institute,  University  of  Bergen, 
Bergen,  Norway,  Tj0tta  was  a  Visiting  Research  Fellow  at  ARL:UT  during  the 
period  1987-89.  Although  only  partly  supported  by  ONR  funds,  he  contributed 
significantly  to  the  projects  of  Cotaras  and  Ten  Cate. 

7.  Wayne  M.  Wright  —  Wright  (Physics  Department,  Kalamazoo  College,  Kala¬ 
mazoo,  Michigan)  consulted  annually  for  NAD/ARL:UT  during  1984-94.  He 
also  spent  his  1989-90  sabbatical  leave  (9  months)  at  NAD/ARL:UT  as  a  Vis¬ 
iting  Research  Fellow.  Although  this  appointment  required  no  salary  support 
from  ONR  funds,  Wright  contributed  significantly  to  ONR  supported  projects 
during  the  period. 
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3.  PROJECTS 


For  purposes  of  this  report,  the  broad  headings  below  are  used  to  classify  the 
research  that  was  done.  The  categories  are  not  mutually  exclusive,  however,  and  the 
classification  for  some  projects  is  arbitrary. 


1.  Nonlinear  Ocean  Acoustics 

2.  Propagation  Problems 

3.  Reflection,  Refraction,  Interaction,  and  Scattering 

4.  Problems  with  Application  to  Medical  Ultrasonics 

5.  Laboratory  Equipment 

6.  Miscellaneous  and  General 

3.1  Nonlinear  Ocean  Acoustics 


For  background  see  the  final  report  for  Contract  N00014-75-C-0878  (86-1,  pp.  42- 
44).  Two  projects  begun  during  the  earlier  period  were  completed.  First  was  a 
study  of  the  effect  of  nonlinear  propagation  distortion  on  pulses  that  travel  great 
distances  in  a  stratified  ocean  (84-4,  85-2,  85-3,  85-7,  86-2,  86-3,  87-1).  The  work 
was  carried  out  primarily  by  Cotaras,  but  Morfey  (and  to  a  minor  extent  R.  Buckley) 
also  contributed.  The  second  project  had  to  do  with  ocean  properties  important 
in  nonlinear  ocean  acoustics,  principally  the  coefficient  of  nonlinearity  /5.  Empirical 
relations  were  developed  that  show  the  dependence  of  /3  and  a  related  coefficient  on 
pressure,  temperature,  and  composition,  for  both  fresh  water  and  seawater  (91-8, 
93-4). 


3.2  Propagation  Problems 


A  wide  variety  of  projects  fall  in  this  category.  They  are  listed  here  roughly  in 
chronological  order. 
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a.  Finite-amplitude  propagation  in  bulk,  air-filled  porous  material.  For  a  summary 
of  previous  work  on  this  project,  see  86-1,  pp.  44-47,  and  the  journal  article  87- 
3.  The  project  was  completed  by  Nelson  early  in  the  period  covered  by  this 
report.  Nelson’s  technical  report  (85-1)  covers  his  master’s  thesis  (84-8)  but 
also  contains  a  chapter  of  new  material.  See  also  85-4  and  90-6. 

b.  Nonlinear  propagation  effects  in  a  rectangular  waveguide.  The  principals  here 
were  Ten  Cate  and  Hamilton.  The  work  was  an  outgrowth  of  Hamilton’s 
Ph.D.  research  (84-5,  84-6)  and  Ten  Cate’s  M.S.  research  (84-2).  Collaboration 
between  the  two,  begun  in  1982—83  while  both  were  graduate  students,  resumed 
when  Hamilton  rejoined  NAD/ARL:UT  as  a  postdoctoral  fellow  (1984-85)  and 
Ten  Cate  returned  in  1985  to  begin  his  doctoral  studies.  The  work  was  both 
theoretical  and  experimental.  Two  different  problems  were  investigated: 

•  Excitation  of  the  waveguide  in  the  1,0  mode  by  a  single-frequency  source, 
which  produces  a  family  of  harmonic  distortion  components  (85-5,  86-4, 
88-6). 

•  Interaction  of  two  noncollinear  waves,  one  of  low  frequency  in  the  0,0  mode, 
the  other  of  high  frequency  in  the  1,0  mode  (85-5,  85-9,  86-4,  87-2). 

c.  Intensity  and  absorption  of  finite-amplitude  waves.  This  work  by  Blackstock 
was  stimulated  by  the  need  to  predict  an  effective  measure  of  absorption  for 
finite-amplitude  ultrasound  in  tissue  (88-3,  90-5). 

d.  Waves  in  a  periodic  medium.  Bradley’s  pioneering  work  on  this  topic  grew  out 
of  his  desire  to  investigate  the  interaction  of  nonlinear  propagation  distortion 
with  dispersion.  However,  the  vehicle  he  chose  for  the  study,  propagation  in 
a  periodic  waveguide,  led  him  into  a  wide-ranging  exploration  of  small-signal 
acoustical  propagation  in  a  periodic  medium,  including  general  aspects  of  phase 
and  group  velocity,  as  well  as  nonlinear  propagation  effects.  The  work  was  both 
theoretical  and  experimental.  It  resulted  in  several  oral  presentations  (89-4, 
90-3,  91-2,  93-8),  a  conference  proceedings  paper  (90-7),  Bradley’s  master’s 
thesis  (90-11,  91-4)  and  doctoral  dissertation  (93-1),  and  three  journal  articles 
(94-3,  94-4,  95-6).  Largely  because  of  the  quality  and  success  of  his  graduate 
work,  Bradley  was  awarded  the  Acoustical  Society’s  F.  V.  Hunt  Postdoctoral 
Fellowship  for  1993-94. 

e.  Momentum  equation  for  plane  progressive  waves.  This  was  a  brief  study  carried 
out  by  Hamilton  and  Blackstock  to  publicize  a  little-known  fact  about  nonlin¬ 
ear  propagation  distortion:  for  plane  progressive  waves  nonlinear  behavior  is 
associated  with  nonlinear  terms  in  the  continuity  equation  and  the  equation  of 
state;  second-order  terms  disappear  from  the  momentum  equation  (90-9). 

f.  Finite-amplitude  radiation  from  a  baffled  piston.  As  part  of  his  project  on  scat¬ 
tering  of  sound  by  sound  (see  Sec.  3.3.c  below),  Ten  Cate  made  a  very  careful 
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and  thorough  study  of  the  field  of  a  single  piston  source  driven  hard  at  a  single 
frequency  (91-7,  93-6).  Although  good  measurements  of  finite-amplitude  radi¬ 
ation  from  a  piston  have  been  reported  by  previous  investigators,  Ten  Cate’s 
results  are  in  a  class  by  themselves.  Dynamic  range  of  his  measurements  fre¬ 
quently  exceeds  100  dB,  and  comparison  between  measurement  and  theory  (the 
“Bergen  code”  is  shown  in  very  fine  detail. 

g.  Finite-amplitude  waves  in  a  medium  having  a  distribution  of  relaxation  pro¬ 
cesses.  Li’s  project  was  motivated  by  concern  about  biomedical  use  of  high- 
intensity  ultrasound.  (Despite  the  obvious  connection  with  Sec.  3.4,  the  re¬ 
search  is  listed  here  because  of  its  fundamental  nature  and  also  because  of  its 
possible  application  to  other  media,  such  as  sediments.)  The  study  was  twofold. 
First,  Li  developed  an  equation  to  model  finite-amplitude  propagation  in  tissue. 
The  absorption  coefficient  of  tissue  is  known  empirically  to  vary  with  (angular) 
frequency  as  w”,  where  n  is  generally  in  the  range  1.0-1. 2.  Carstensen  explained 
the  peculiar  frequency  dependence  by  modeling  tissue  as  a  medium  having  a 
broad  distribution  of  relaxation  processes.^  Using  Carstensen’s  model,  Li  de¬ 
rived  a  nonlinear,  integro-differential  equation  (a  generalized  Burgers  equation) 
for  propagation  of  finite-amplitude  waves  in  a  thermoviscous,  multiply  relaxing 
fluid.  Some  asymptotic  solutions  and  some  numerical  solutions  were  found  (92- 
4,  93-7,  93-9).  In  the  second  part  of  the  study,  Li  applied  his  generalized  Burgers 
equation  to  the  problem  of  enhanced  absorption,  and  concomitant  temperature 
rise,  in  tissue  exposed  to  intense  ultrasound  (93-7,  93-9,  94-1).  He  used  the 
bioheat  equation  to  calculate  the  increase  in  temperature. 

h.  Nonlinear  distortion  in  propagation  of  aircraft  noise.  This  was  a  project  begun 
with  ONR  support  but  soon  shifted  to  NASA  support  (95-2). 


3.3  Reflection,  Refraction,  Interaction,  and  Scattering 

a.  Angular  dependence  of  the  coefficient  of  nonlinearity  p.  For  ordinary  propa¬ 
gation  0  depends  only  on  the  material  property  B/A  of  the  fluid  (the  formula 
is  =  1  4-  BI2A).  When  two  noncollinear  wave  fields  interact,  however,  the 
nonlinear  effect  of  one  field  on  the  other  depends  on  the  angle  9  at  which  they 
intersect.  Our  work  on  this  problem  began  under  the  previous  contract,  N00014- 
75-C-0867;  see  84-2  and  also  86-1,  pp.  24-29.  The  topic  is  also  closely  related  to 
the  waveguide  work  described  in  Sec.  3.2.b  above;  see  in  particular  85-9.  Hamil¬ 
ton  and  Blackstock  wrote  a  journal  article  (88-1)  to  discuss  the  theoretical  basis 
and  validity  of  the  expression  /3  =  cos0  -f-  B/2A,  which  often  appears  in  the 
literature. 

b.  Reflection  and  refraction  at  a  plane  interface  between  two  fluids.  This  topic 
was  Cotaras’s  doctoral  project.  In  the  beginning  the  goal  was  (1)  to  rederive 
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c. 


Snell’s  law  and  the  law  of  specular  reflection  for  plane  finite-amplitude  waves 
that  are  obliquely  incident  on  a  plane  interface  between  two  fluids,  and  (2)  to 
test  the  resulting  predictions  by  experiment.  Initially  an  apparently  successful 
derivation  of  a  new  form  of  Snell’s  law  was  obtained  by  physical  arguments  and 
various  approximations,  and  a  test  experiment  was  planned  (see  87-4,  pp.  3-7, 
and  also  87-7).  However,  the  derivation  was  found  to  be  too  simple  (and  the  ap¬ 
proximations  questionable),  and  very  substantial  difficulties  were  encountered 
in  the  experiment  (88-5,  pp.  2-4).  The  task  was  therefore  changed  to  a  theo¬ 
retical  investigation  of  oblique-angle  refraction/reflection  for  weakly  nonlinear 
waves,  but  with  all  sources  of  nonlinear  effects  taken  into  account,  including 
motion  of  the  interface,  exact  source  conditions,  and  nonlinear  relations  among 
various  field  variables.  The  results  are  exceedingly  complicated.  However,  it 
appears  that  SneU’s  law  and  specular  reflection  continue  to  hold,  at  least  to 
second  order,  for  finite-amplitude  waves  (89-1,  89-7).  The  Tjpttas  were  very 
helpful  to  Cotaras  on  theoretical  aspects  of  this  project. 

Scattering  of  sound  by  sound.  Ten  Cate’s  doctoral  research  was  intended  to  test 
new  predictions  by  the  Tjpttas^  and  Hamilton  and  Darvennes^  on  scattering  of 
sound  produced  when  two  noncollinear  beams  of  different  frequency  interact. 
In  1956  Ingard  and  Pridmore-Brown  claimed®  to  have  observed  scattered  sound 
outside  the  region  where  the  two  beams  interact.  The  following  year,  and  many 
times  since,  Westervelt,®  on  the  basis  of  theoretical  considerations,  categori¬ 
cally  denied  the  existence  of  scattered  sound  outside  the  interaction  region. 
Ever  since,  the  question  has  been  highly  controversial.  Ten  Cate  performed 
new  experiments,  which  were  preceded  and  complemented  by  new  theoretical 
work  with  the  Tjpttas  (89-8,  89-9,  91-1).  Two  sets  of  experiments  were  done. 
The  first  was  on  “self-scattering”  of  sound  by  sound,  which  occurs  when  only 
one  source,  driven  at  a  single  frequency,  is  used.  The  self-scattering  manifests 
itself  as  extra  sidelobes,  called  fingers,  in  the  beam  patterns  of  the  second- 
harmonic  and  higher  distortion  components.  Since  the  fingers  exist  outside  the 
interaction  region  (in  this  case  the  beam  pattern  of  the  fundamental  defines  the 
interaction  region),  they  qualify  as  scattered  sound.  Ten  Cate’s  measurements 
(Sec.  3.2. f  above)  were  outstandingly  successful  (91-7,  92-1,  93-6).  No  one  before 
had  so  carefully  measured  the  beam  patterns  and  compared  them  with  theory, 
particularly  beyond  the  second  harmonic.  Moreover,  he  was  the  first  to  mea¬ 
sure  propagation  curves,  both  for  fingers  and  the  product  pattern  lobes,  and  to 
quantify  their  different  diminutions  with  distance.  The  second  set  of  measure¬ 
ments  was  on  crossed-beam  scattering  of  sound  by  sound.  Unfortunately,  they 
were  not  so  successful.  Ideally,  to  avoid  the  presence  of  primary  beam  sidelobes, 
which  interfere  with  observation  of  scattered  sound,  the  experiment  should  be 
done  with  two  Gaussian  sources.  Only  one  approximately  Gaussian  source  was 
available,  and  it  did  not  work  well.  The  final  results  were  inconclusive  (92-1, 
94-9).  Any  future  attempt  to  observe  crossed-beam  scattering  should  be  done 
with  two  true  Gaussian  sources. 
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d.  Finite-amplitude  waves  in  a  three-layer  complex.  Yazdi,  who  was  interested  in 
biomedical  ultrasonics,  was  given  this  project  because  of  its  potential  application 
to  tissue.  The  project  was  designed  to  use  finite-amplitude  sound  in  another 
traditional  problem  drawn  from  classical  linear  acoustics:  transmission  of  sound 
from  medium  1  into  medium  2  and  thence  into  medium  3.  Frequently  media  1 
and  3  are  the  same,  for  example,  when  the  application  is  to  transmission  loss 
through  a  panel  or  wall.  Although  several  aspects  of  this  problem  could  be 
explored,  the  one  selected  for  Yazdi  was  a  proposal  (and  theoretical  prediction) 
by  Nazarov^  for  measuring  the  coefficient  of  nonlinearity  /3.  In  this  arrangement 
a  layer  of  material  (medium  2)  in  a  water  bath  (media  1  and  3)  is  exposed 
to  normally  incident  plane  waves  of  finite  amplitude.  The  source  frequency  is 
chosen  so  that  the  layer  of  material  is  transparent  to  the  incident  sound,  which  is 
composed  of  the  fundamental  and  also  the  second-harmonic  distortion  generated 
in  the  water.  In  this  case  the  only  second-harmonic  signal  reradiated  back  to  the 
source,  which  also  serves  as  a  receiver,  is  that  from  the  standing  wave  second- 
harmonic  field  generated  in  the  layer  of  material.  Measurement  of  the  reradiated 
signal  allows  0  for  the  material  to  be  determined  (93-2,  pp.  16-19).  For  initial 
simplicity  a  steel  plate  was  used  for  the  layer  of  material.  Unfortunately  before 
many  measurements  could  be  taken,  the  experiment  was  stopped  because  Yazdi 
transferred  to  a  more  traditional  biomedical  engineering  program  within  the 
University.  He  had,  however,  made  a  very  valuable  contribution  by  designing 
and  building  a  water  tank  and  computer  controlled  positioning  apparatus  (see 
Sec.  3.5.b  below).  Although  intended  for  Yazdi’s  experiment,  the  facility  has 
proved  extremely  useful  for  several  other  NAD/ARL:UT  projects. 

e.  Self-refraction  in  the  field  of  a  paraboloidal  reflector.  This  was  Gelin’s  project. 
Self-refraction  is  bending  of  rays  due  solely  to  finite-amplitude  effects,  as  op¬ 
posed  to  ordinary  refraction,  which  is  caused  by  sound  speed  variation  in  an 
inhomogeneous  fluid.  The  propagation  speed  of  any  point  on  a  finite-amplitude 
wave  is  ^  =  Cq  +  where  cq  is  the  small-signal  sound  speed  and  u  is  the 
particle  velocity  at  the  point.  Because  of  the  dependence  of  ^  on  u,  any  vari¬ 
ation  of  the  particle  velocity  along  a  wavefront  causes  a  change  in  direction  of 
the  wavefront,  i.e.,  refraction.  Self-refraction,  first  proposed  by  Whitham  for 
shock  waves,®’®  is  difficult  to  isolate  in  acoustical  experiments  because  it  is  often 
mixed  in  with  other  effects,  for  example  diffraction  in  the  radiation  field  of  a 
piston.  However,  the  paraboloidal  reflector,  with  a  spark  source  at  the  focus, 
offers  an  ideal  means  for  observing  self-refraction:  the  (linear  theory)  plane 
wave  reflected  field  is  inhomogeneous,  that  is,  the  amplitude  varies  along  the 
wavefronts.  As  a  result,  the  wavefronts  are  predicted  to  become  curved,  concave 
for  negative  parts  of  the  waveform,  convex  for  positive  parts.  Gelin  carried  out 
the  experiment  in  air.  He  made  a  large  number  of  microphone  measurements 
over  the  entire  reflected  field.  Some  peculiar  effects  were  observed  near  the  axis, 
probably  due  to  the  effect  of  the  spark  source  (either  local  heating  of  the  air  or 
diffraction  by  the  electrodes),  which  affected  the  paraxial  reflected  rays.  Despite 


11 


these  difficulties,  self-refraction  was  very  clearly  observed  (94-6,  94-7,  95-1). 


3.4  Problems  with  Application  to  Medical  Ultrasonics 

a.  Lithotripsy.  Lithotripsy  is  the  use  of  focused  shock  waves  to  break  up  kidney 
stones.^"  In  the  most  common  lithotripter,  the  Domier  machine,  the  shock  wave 
source  is  an  ellipsoidal  reflector  with  an  electric  spark  at  its  near  focus.  Both 
source  and  patient  are  in  a  water  bath.  The  patient  is  positioned  so  that  the 
kidney  stone  is  located  at  the  far  focus  of  the  ellipsoid.  The  pressure  waveform 
at  the  far  focus  is  a  very  short  positive  spike,  of  amplitude  as  high  as  100  MPa 
and  duration  about  1  (jls,  followed  by  a  negative  tail  of  much  longer  duration  and 
magnitude  of  order  10  MPa.  Several  elements  -  very  high  amplitude,  focusing, 
diffraction,  and  cavitation  -  make  lithotripsy  a  fascinating  problem  in  nonlinear 
acoustics.  The  physics  of  lithotripsy  prompted  our  general  interest,  beginning 
in  1987  (87-6,  88-2,  88-4).  Work  on  specific  topics  followed;  see  Secs.  3.4.b  and 
3.4.d  below. 

b.  Ellipsoidal  focusing  in  air.  This  task  was  set  up  to  study  the  interaction  of  finite- 
amplitude  propagation,  focusing,  and  diffraction  in  the  field  of  an  ellipsoidal 
reflector.  Although  lithotripsy  was  the  intended  application,  the  measurements 
were  made  in  air  because  of  our  previous  experience  with  airborne  pulses  gener¬ 
ated  by  sparks  (86-1,  pp.  29-39).  S.  T.  W.  Cheng  (M.E.  undergraduate),  Wright, 
and  Blackstock  carried  out  the  study.  Four  different  ellipsoidal  reflectors  were 
machined.  Just  as  in  lithotripsy,  a  spark  located  at  the  near  focus  generated 
the  shock  pulse.  Although  the  waveform  in  air,  an  N  wave,  is  different  from 
that  in  water  (a  sharp  positive  peak  followed  by  a  long  shallow  negative  tail), 
the  physics  of  the  reflection,  focusing,  and  diffraction  processes  is  the  same. 
Pressure  signatures  of  the  reflected  wave  were  measured  at  a  sequence  of  points 
along  the  axis:  (1)  in  the  prefocus  region,  where  the  waveform  observed  was 
still  an  N,  (2)  in  the  focal  region,  where  the  waveform  was  U-shaped,  and  (3)  in 
the  postfocus  region,  where  the  waveform  was  a  “sagging  N”  (89-6,  pp.  9-11). 
In  linear  theory,  passage  through  a  three-dimensional  focus  is  accompanied  by 
waveform  inversion,  that  is,  a  180°  phase  shift.  It  is  clear,  therefore,  that  strong 
nonlinear  distortion  in  the  focal  region  and  beyond  greatly  modified  the  phase 
change  usually  associated  with  focusing.  Other  measurements  and  observations 
were  made  (89-2,  89-3,  91-6). 

c.  Negative  pulse  generation  by  diffraction  and  other  means.  This  was  Bailey’s 
M.S.  research  project.  An  isolated  negative  pressure  pulse  in  water  was  needed 
for  a  medical  ultrasonics  experiment.  Bailey  began  with  an  electric  spark,  which 
produces  an  almost  pure  positive  pulse  (the  long  negative  tail  that  follows  is 
frequently  so  weak  as  to  be  lost  in  the  noise).  Next  he  inverted  the  positive  pulse 
by  one  of  two  methods:  (1)  reflection  from  the  water  surface,  or  (2)  diffraction 
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from  a  circular  aperture  (the  edge  wave  on  axis  behind  the  aperture  is  an  in¬ 
verted  replica  of  the  incident  wave).  Without  further  steps,  the  biomedical 
target  would  then  have  been  exposed  to  the  positive  direct  pulse  as  well  as  the 
desired  negative  pulse.  To  block  out  the  direct  pulse,  Bailey  made  use  of  the  fact 
that  scattering  from  an  object  having  a  ragged  edge  produces  an  incoherent, 
very  weak  diffracted  signal.  For  method  (1),  reflection  from  the  water  surface, 
he  interposed  a  ragged-edge  corprene  barrier  (either  a  plate  or  a  disk)  between 
the  spark  source  and  the  target.  For  method  (2),  the  aperture,  he  placed  a  small 
ragged-edge  corprene  disk  in  the  aperture.  In  both  cases  the  blocking  object  got 
rid  of  the  direct  wave  and  yet,  because  of  the  ragged  edge,  introduced  no  new 
coherent  edge  wave.  The  work  was  described  in  three  oral  reports  (92-3,  94-5, 
95-7),  one  conference  proceedings  (93-3),  and  Bailey’s  thesis  (94-2).  Method 
(1),  with  a  plate  used  to  block  out  the  direct  wave,  was  used  in  a  biomedical 
experiment  conducted  at  the  University  of  Rochester.  Mouse  lungs  and  fruitfly 
larvae  were  exposed  to  isolated  positive  pulses  and  to  isolated  negative  pulses 
(94-8,  96-1).  The  positive  pulses  were  found  to  be  at  least  as  damaging  as  the 
negative  pulses.  This  is  an  important  finding.  Heretofore  investigators  had 
assumed  that  tissue  damage  is  due  primarily  to  negative  pressure. 

d.  Intensified  cavitation  produced  by  a  duo  ellipsoidal  reflector  system.  Bailey’s 
doctoral  research  project  (at  this  writing  about  half  finished)  is  built  around 
a  two-ellipsoid  system  that  is  designed  to  enhance  the  violence  of  cavitation 
collapse.  The  work  is  most  concisely  described  by  the  following  abstract  of  a 
paper  to  be  presented  at  an  upcoming  meeting  (96-2): 

“An  underwater  bubble  is  well  known  to  grow  in  response  to  a  strong  nega¬ 
tive  acoustic  pulse  and  then  collapse  because  of  inertial  forces.  Here  we  show 
that  adding  an  auxiliary  positive  pulse,  after  collapse  begins,  intensifies  the 
collapse.  The  negative-then-positive  pulse  sequence  is  produced  by  two  ellip¬ 
soidal  reflectors,  each  with  an  electrical  spark  at  its  near  focus  /i  and  beamed 
so  that  they  share  a  common  second  focus  f2-  The  negative  pulse  is  produced 
by  a  polyurethane  (pressure  release)  ellipsoid,  the  positive  by  a  brass  (rigid) 
ellipsoid.  A  timing  circuit  controls  the  delay  between  the  pulses.  Cavitation  is 
recorded  by  pitting  (caused  by  bubble  collapse)  of  an  aluminum  foil  membrane, 
which  is  centered  at  /a  and  lies  coplanar  with  the  two  crossed  beams.  When  the 
brass  reflector  is  fired  alone,  a  narrow  path  of  ~l-mm  diameter  pits  appears. 
Firing  the  polyurethane  reflector  alone  yields  more  widespread,  ~0.Tmm  di¬ 
ameter  pits.  When  both  are  fired,  a  pitted  X  pattern  shows  the  position  of  the 
two  beams.  If  the  delay  between  the  two  firings  is  2-6  fj,s,  the  intersection  of 
and  the  centerline  between  the  paths  erupts  with  deep  pits.  [Work  supported 
by  ONR.]” 
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3.5 


Laboratory  Equipment 


a.  Nonlinear  acoustics  laboratory  upgrade.  A  substantial  upgrade  of  the  facili¬ 
ties  for  nonlinear  acoustics  experiments  took  place  during  the  two-year  period 
1988-90.  Support  for  the  project  came  jointly  from  ONR  and  ARL:UT  IR&D. 
Cotaras  and  Bradley  were  largely  responsible  for  the  design,  equipment  acqui¬ 
sition,  and  integration  of  a  new  measurement  system,  which  made  possible  a 
Avide  variety  of  computer  controlled  experiments.  See  89-6  and  90-4  for  a  de¬ 
tailed  description.  The  new  system  played  an  extremely  important  role  in  all 
subsequent  experiments,  e.g.,  those  described  in  Secs.  3.2.d,  3.3.d,  3.3.e,  3.4.b, 
3.4.C,  and  3.4.d,  and  in  sonic  boom  model  experiments  supported  by  NASA.'^ 

b.  Ultrasonics  tank  with  computer  controlled  positioner.  In  1990-91  Yazdi  greatly 
expanded  the  capability  of  the  measurement  system  described  above  by  adding 
a  water  tank  and  precision  positioning  system,  which  he  designed  and  built 
himself.  NIH  provided  funds  for  the  motor  controller/driver  and  the  motors, 
which  were  the  most  expensive  components.  The  tank  and  positioner,  or  in  some 
cases  the  positioner  by  itself,  were  heavily  used  in  several  subsequent  projects, 
for  example,  those  described  in  Secs.  3.3.d,  3.3. e,  3.4.c,  and  3.4. d,  and  the  sonic 
boom  experiments.^^ 


3.6  Miscellaneous  and  General 

a.  Subharmonic  generation  and  chaos  in  an  acoustical  resonance  tube.  This  started 
out  as  Ten  Cate’s  doctoral  project.  The  premise  was  that  a  closed  end  resonance 
tube  driven  in  one  of  its  higher  plane  wave  modes  would  excite  subharmonics 
(in  the  lower  modes  of  the  tube).  If  the  tube  were  driven  hard  enough,  chaos 
would  result.  After  careful  study,  however.  Ten  Cate  found  that  the  resonance 
tube  is  not  a  good  candidate  for  chaotic  response.  A  major  problem  is  that 
raising  the  drive  amplitude  rapidly  increases  the  losses  and  thus  lowers  the  Q 
of  the  system.  Since  the  resonance  response  curve  cannot  bend  over  far  enough 
to  produce  hysteresis,  bifurcation  cannot  occur.  Ten  Cate  wrote  an  excellent 
white  paper  on  the  subject;  see  86-6,  Appendix  A. 

b.  Thermoacoustics.  Wright’s  work  on  the  acoustical  response  of  a  cylindrical 
cavity  to  a  modulated  heat  source  (a  current  carrying  wire)  was  reported  in  a 
journal  article  (87-5). 

c.  Reflection  and  transmission  of  a  spherical  wave  at  a  concentric  spherical  in¬ 
terface.  Reflection  and  transmission  of  plane  waves  normally  incident  on  a 
plane  interface  between  two  lossless  fluids  (impedance  Z\  =  p\C\  on  the  inci¬ 
dent  wave  side,  Z%  —  P2C2  on  the  transmitted  wave  side)  is  a  classical  subject. 
Blackstock  and  Morfey  investigated  the  analogous  problem  for  spherical  waves 
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(91-3).  Among  the  results:  (1)  The  reflection  and  transmission  coefficients  R 
and  T  are  complex  and  frequency  dependent;  the  coefficients  for  a  convex  in¬ 
terface  are  the  complex  conjugates  of  those  for  a  concave  interface.  (2)  At  high 
frequency  R  and  T  approach  their  plane  wave  values;  at  very  low  frequency 
J?  — >•  —  1  (unless  Pi  =  P2).  (3)  Perfect  transmission  requires  both  pi  =  p2  and 
Cl  =  C2,  not  just  ^1  =  ^2- 

d.  Acoustitron.  The  acoustitron  is  a  toroidal  (doughnut  shaped)  waveguide  driven 
from  two  or  more  (side-mounted)  sources  that  are  phased  to  assure  unidirec¬ 
tional  propagation  around  the  waveguide.  When  the  frequency  is  such  that 
the  (axial)  circumference  of  the  toroid  is  an  integral  number  of  wavelengths, 
resonance  occurs.  Resonance  in  a  progressive  wave  fleld  is  very  unusual  and 
is  what  makes  the  acoustitron  interesting.  The  high  amplitude  that  may  be 
obtained  also  makes  it  an  attractive  candidate  for  nonlinear  acoustics.  Bailey 
briefly  investigated  the  acoustitron  as  a  possible  doctoral  project.  Two  theoret¬ 
ical  results  were  obtained:  (1)  an  expression  for  the  gain  that  can  be  achieved 
from  the  progressive  wave  resonance,  (2)  a  calculation  of  the  distorted  wave¬ 
form  that  develops  because  of  nonlinear  propagation  effects.  Some  preliminary 
experiments  were  also  done.  See  95-2,  pp.  8-11,  for  details.  In  the  end  Bailey 
preferred  the  duo  ellipsoid  project  (Sec.  3.4. d  above)  for  his  doctoral  work,  and 
the  acoustitron  project  was  set  aside  for  a  future  student. 

e.  General.  Other  work  associated  in  some  way  with  Grant  N00014-89-J-1109, 
Contract  N00014-84-K-0574,  or  one  of  its  predecessors  is  as  follows:  several  re¬ 
views  of  nonlinear  acoustics  (86-5,  89-5,  90-2,  90-10,  95-4,  95-5),  journal  articles 
on  earlier  work,  either  appearing  for  the  first  time  or  reprinted  in  books  (84-9, 
84-10,  85-6),  historical  sketches  (84-7,  92-2),  a  paper  on  acoustics  education 
(88-7),  and  a  book  (90-8). 
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4.  SUMMARY 


Research  on  nonlinear  acoustics  was  performed  during  the  period  1984-95.  The 
major  areas  included  in  the  report  are  as  follows:  nonlinear  ocean  acoustics;  propaga¬ 
tion  problems;  reflection,  refraction,  interaction,  and  scattering;  problems  associated 
with  medical  ultrasonics;  and  miscellaneous  topics.  The  research  led  to  flve  mas¬ 
ter’s  degrees  and  four  doctorates  among  seven  students.  One  other  doctorate  was  in 
progress  at  the  time  the  grant  ended.  Also  participating  in  the  research  were  one 
postdoctoral  fellow,  four  visiting  faculty  members  from  other  universities,  and  a  vis¬ 
itor  from  a  Canadian  naval  laboratory.  Among  the  public  disclosures  of  the  research 
were  15  journal  articles  (and  two  more  submitted  but  not  yet  published) ,  10  papers 
in  conference  proceedings,  6  book  chapters,  39  papers  presented  orally  at  national 
and  international  meetings,  9  theses  and  dissertations,  and  18  technical  reports. 
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